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Abstract - It is necessary to examine the contact stresses 
developed between the contacting bodies. The contact stress 
depends on the material property of the contact bodies, 
clearance between them and load acting on it. The combination 
of contacting materials and contact conditions are critical. 
The effect of clearance and load on the contact conditions of 
cylindrical conformal contact is examined for same and varies 
values of Young's modulus of the two contacting bodies. Finite 
element analysis is used to study cylindrical conformal 
contacts with the models and is validated by comparison with 
known analytical solutions. These results can be applied to 
any problem of conformal contact in nature such as actuators, 
hip-joint, piston barrel, cylinder block and the piston ball, 
slipper seat of a water-lubricated axial piston pump and for 
solid structural bearing with shaft. 

Index Terms - Contact Stress, Conformal Cylindrical Contact, 
Contact Angle, Contact Clearance, Half Space. 

I. Introduction 

The determination of pressure at a contact in a machine 
part is important because contact stresses frequently lead to 
failure by seizure, wear or fatigue. Due to a small radial 
clearance and a heavy load application in certain mechanical 
joints a conformal (Non-Hertizian) type of contact exists. 
Under the application of the load the size of the contact area 
grows rapidly and may become comparable with the significant 
dimensions of the contacting bodies. A pin in a hole with a 
small clearance is an example. When the area of contact 
occupies an appreciable fraction of the circumference of the 
hole neither the pin nor the hole can be regarded as an elastic 
half-space so that the Hertz treatment is invalid. But still it 
can be considered as half-space for the purpose of calculating 
elastic deformations and stresses. 

Many problems of contact mechanics are solved using 
the complex integration are given by N.I.Muskhelishvili in 
the year 1963[1]. The finite element analysis and approximate 
model for cylindrical joints with clearances are given by Cai- 
Shan Liu, Ke Zhang and Rei Yang in 2006 [2] where By 
introducing some appropriate assumptions and analyzing the 
FEM numerical results, such as the contact area, the pressure 
distribution, and the maximum sustainable load, an 
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approximate model for the contact problem of cylindrical joints 
with clearances is developed through modeling the pin as a 
rigid wedge and the elastic plate as a simple Winkler elastic 
foundation [3]. 

There are numerous methods to solve contact problems. 
Some of the different methods of contact stress evolution 
are analytical closed form solution, numerical solutions and 
experimental methods. It is highly desirable to obtain complete 
closed-form analytical solutions to problems in contact 
mechanics. An analytic result means, that the full internal 
stress-displacement field for the given situation is expressed 
in a closed form in terms of known elementary (polynomial, 
exponential, etc.) and special (elliptical, Bessel, etc.) functions 
[4] . The numerical techniques are powerful tools to analyze 
the contact problem because of their flexibility and ability to 
model all complications involved with the analysis of such 
problems. 

Recently, several contact algorithms have been proposed 
and incorporated into commercially available Finite Element 
software. In the contact zones existing in various mechanisms, 
stresses due to quasi-static loads can lead to friction and to 
high stresses thereby inducing fatigue and erosion of 
surfaces. In order to design these elements it is necessary to 
predict contact actions and stresses. Exact solutions of 
contact problems are available for ideal conditions [5,6,7]. 

The combined effects on the stress concentrations of the 
shaft radius, the interference, and Young's modulus of shaft 
and hub with the aid of finite elements; various design plots 
are compiled that report the elastic stress concentrations, 
contact pressure within the hub versus the normalizing 
parameter like contact angle and clearance between the shaft 
and hub. In most of the actual situations analytical solutions 
don't exists and approximate solutions have to be obtained, 
using either Finite Element method formulation or a variation 
formulation. 

II. THEORETICALFORMULATION 

Steumann found the distribution of pressure P (x) for 
profiles having the form A x 2n . For the two dimensional half- 
space problems in which the displacements are specified over 
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the interval -a < x < a , the singular integral equation is 
obtained as 

m 

Here F(s) is the component of traction. The resulting equation 

is 
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Here the Hertz theory correspond to n= 1 . For higher values 
of n the pressure has its maximum values away from the centre 
of the contact. Persson assumed that the contact surface 
was cylindrical and formulated the integro-differential 
equation and determined the analytical contact pressure 
distribution for the plane stress case. The contact stress 
distribution for the plane stress condition is given by Persson 
by the expression, 



I>(t) = 



2F<c 2 -q 2 )^ 



ttR<l +q 3 )(c 2 



2ieRc : 



xln 



(c 2 +l)^ + (c 2 -q 2 )^ 



<5) 



(c 2 +l)^-(c 2 -q 2 )^ 
The expression for the normalized maximum pressure at the 
centre of contact is calculated from the eq-5. When (j) = , 
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An expression relating the contact angle a with the load F, 
radial difference CL, Young's modulus E is 
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The eq-7 is solved by Complex integration method and the 
resulting equation is given by 
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Where c 2 = x , c = tanl j j > q = tan| -^ | ; and 

F = Applied load, 

^ =The angular position, 

1(X =The angle of contact and 
R = contacting radius. 
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Figure 1. Cylindrical Conformal Contact 



m. FINITE ELEMENT FORMULATION 

Considering the symmetry of the problem the quarter of a 
cylinder inside a semi-infinite cylindrical cavity is taken and 
the non-contacting half has been neglected as shown in 
figure-2. Where R = 10mm, CL= 1, 0.5, 0.25 mm, X= 200 mm, 
Y = 200 mm. The dimension of the cavity is taken as large to 
satisfy the half space condition. 
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Figure 2. Geometry of the Model and mesh pattern 

The bottom edge of the cavity is fixed in all direction and 
the left side edge is given symmetrical boundary condition. 
The plane stress analysis is done using Finite Element 
Analysis Software. The mesh near the contact area is refined 
sufficiently and the aspect ratio is maintained nearly equal to 
ideal. The analysis is carried out considering same material 
for the cylinder and cavity by taking the young's modulus of 
2.1 x 10 5 N/mm 2 and Poisson's ratio of 0.3, and different material 
for the cylinder and cavity. The modulus index is defined for 
the different material for the cylinder and the cavity as the 
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ratio of the Young's modulus of the cylinder and that of the 
cavity and is denoted by 'n'. 



E,*E 2 



(9) 



In modeling the contact problems in the static analysis 
rigid body motion is the common problem. So to remove this 
problem modeling can be done in such a way that the contact 
bodies are just touching at a single point. The numbers of 
nodes near the contact area that may probably be come in 
contact is identified on two bodies and has been created the 
contact pair using those nodes. Both the bodies are 
considered as flexible one. The two node two dimensional 
surface-to-surface contact element is used with target element 
to create contact pairs. These contact elements are created 
on the exiting structural elements and the two dimensional 
surface-to-surface contact definitions are used for the 
analysis. 

IV. RESULTS AND DISCUSSION 

A. Same Young's Modulus for Two Contacting Bodies 

The Young's modulus of the two contacting bodies is 
taken same. The results are plotted for three different 
clearances varied as 0.25mm, 0.5mm and 1mm at constant 
load. The Contact Pressure along contact surface for the 
various loads is as shown in fig3. 




Figure 3. The Contact Pressure along contact surface for the load 
10KN 





Figure 5. Variation of Contact pressure with the contact angle for 
the load F=10000N and 0.5mm Clearance. 




Figure 6. Variation of Contact pressure with the contact angle for 
the load F=10000N and 1mm Clearance. 

The contact pressure increases with the increase in 
clearance between two contacting bodies and hence 
decreasing the contact angle. More the surface in contact 
more the stress developed due to dynamic friction as shown 
in fig 4 to 6. 
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Figure 7. Von-Mises Stress along the Symmetry Axis of Cavity for 
Different Loads and 0.25mm Clearance. 
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Figure 4. Variation of Contact pressure with the contact angle for 
the load F=10000N and 0.25mm Clearance 
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Figure 8. Von-Mises Stress along the Symmetry Axis of Cavity for 
Different Loads and 0.5mm Clearance 




Figure 9. Von-Mises Stress along the Symmetry Axis of Cavity for 
Different Loads and 1mm Clearance 

The stress varies with the application of load and maxi- 
mum stress is developed at the center of contact region, de- 
creases along symmetric axis as the contact angle increases 
as shown in Fig. 7 to 9. Lesser the clearance more the area of 
contact, so that more the contact angle as the load increases 
and visa-versa. As the load increases the cylinder under- 
goes deformation so that the contact angle increases is shown 
in fig. 10. These plots indicate that as the load increases, 
contact pressure also increases for the contact angle. The 
contact pressure decreases as the contact angle increases. 
This result agrees with 5% error of theory. The Von-Mises 
stress gives the best failure criteria for the steel materials.lt is 
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Figure 10. Variation of Contact Angle with Load for Different 
Clearances 

understood that contact pressure increases with decrease in 
contact angle and as the load increases the Von-Mises stress 
also increases. The contact angle is half the total contact angle 
and is determined theoretically by the equation-8. The varia- 
tion of the contact angle with the load for various values of 
clearances. This shows that as the load increases, the con- 
tact angle increases for the decrease of clearance. 

B. Different Young 's Modulus for Two Contacting Bodies 

Considering different Young's modulus for cylinder and 
cavity (i.e. modulus index), the effect of modulus value on 
the contact pressure is determined. When the young's 
modulus of the cylinder is higher than that of the cavity (i.e. 
n>l) the contact angle increases and the contact stress 
decreases as shown in fig 1 1 . The same trend is observed for 
the increase in load. 

It is noted that as the load increases the contact angle 
increases and contact stress decreases for the case of same 
modulus value by taking the clearance is constant. If one 
decreases the clearance the contact angle increases for the 
constant load case. For a smaller contact angle and contact 
area the maximum pressure is greater for a higher load. All the 
graphs are plotted with well versed computing software 
MATLAB. These results are compromising well with the 
closed form solution obtained. Results well agree with 5% 
error with the theoretical values. So it can be applied to any 
problem of conformal contact in nature considering the 
material properties of contacting bodies, clearance between 
them and the nature of load acting on it. 
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Figure 1 1 . Variation of Contact Pressure and Contact angle for 
different values of Young's Modulus 



CONCLUSION 

It is understood that lesser the clearance between two 
contacting bodies, lesser the contact angle for the constant 
load of application. As the load increases, contact pressure 
also increases for the constant contact angle. The Von-Mises 
stress gives the best failure criteria for the steel materials 
asthe load increases the Von-Mises stress also increases. As 



the load increases the contact angle increases and contact 
stress also increases for the case of same modulus value by 
taking the clearance is constant. If one decreases the clearance 
the contact angle increases for the constant load case. For a 
smaller contact angle and contact area the maximum pressure 
is greater for a higher load. So it is better to have a higher the 
clearance is required to decrease the contact pressure and 
the contact angle is lesser for the increase of the modulus 
ratio. 

REFERENCES 

[1] Muskhelishvili, "Some Basic Problem of the Mathematical 

Theory of Elasticity", N1963. 

[2] Cai-Shan Liu, Ke Zhang and Rei Yang, "The FEM analysis and 

approximate model for cylindrical joints with clearances", Chinese 

Academy of Space Technology, Beijing, China, 2006. 

[3] B. Paul and J. Hashemi, "Contact Pressure on Closely 

Conforming Elastic Bodies", Journal of Applied Mechanics, vol.48, 

pp.543-548, Sept. 1981. 

[4] Hsien H, Chen and Kurt M Marshek, "Effect of Clearance and 

Material Property on Contact Pressure in Two-Dimensional 

Conforming Cylinders", mech. Mach. Theory, Vol.23, No.l, pp.55- 

62, 1988. 

[5] Robert L Jackson and Itzhak Green, "A Finite Element Study 

of Elasto-Plastic Hemispherical Contact against A Rigid Flat", 2005 . 

[6] G. A. Papadopoulos, "Experimental Study of the Load 

Distribution in Bearings by the Method of Caustics and the 

Photoelasticity Method", Journal of Strain Analysis, Vol.40, No.4, 

pp.357-365, 2005. 

[7] M Pau, A Baldi, P. F Orru and F Ginesu, "Experimental 

Investigation on Contact between Cylindrical Conformal Surfaces", 

Journal of Strain Analysis, Vol.39, No.3, pp.315-328, 2004. 

[8] A. Strozzi, A. Baldini, M. Giacopini, E. Bertocchi, L. Bertocchi, 

"Normalization of the stress concentrations at the rounded edges 

of a shaft-hub interference fit", The Journal of Strain Analysis for 

Engineering Design August 2011 vol. 46 no. 6 478-491. 



©2012 AMAE 
DOI:01.DPIE.02.01.508 



26 



^AMAE 



